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Abstract: We consider an Abelian gauge field coupled to a particular truncation of 
Horndeski theory. The Galileon field has translation symmetry and couples non minimally 
both to the metric and the gauge held. When the gauge-scalar coupling is zero the gauge 
held reduces to a standard Maxwell held. By taking into account the symmetries of the 
action, we construct charged black hole solutions. Allowing the scalar held to softly break 
symmetries of spacetime we construct black holes where the scalar held is regular on the 
black hole event horizon. Some of these solutions can be interpreted as the equivalent 
of Reissner-Nordstrom black holes of scalar tensor theories with a non trivial scalar held. 
A self tuning black hole solution found previously is extended to the presence of dyonic 
charge without affecting whatsoever the self tuning of a large positive cosmological constant. 
Finally, for a general shift invariant scalar tensor theory we demonstrate that the scalar 
held Ansatz and method we employ are mathematically compatible with the held equations. 
This opens up the possibility for novel searches of hairy black holes in a far more general 
setting of Horndeski theory. 
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1 Introduction 

The stationary state of a black hole in General Relativity can be determined by a limited 
number of physical charges which are conserved quantities subject to a Gauss law: mass, 
electric and magnetic charge and angular momentum. Quite generically any matter sector, 
which one tries to attach to a black hole is either eaten up and hidden behind the horizon, 
or/and expelled away from the black hole. After relaxation of the system, when we return 
to the stationary state again, the global parameters of the black hole may have changed, 
but ’’hair” does not survive. All properties of a black hole are determined from asymptotic 
measurements. This statement is known as the no-hair paradigm for black holes [1]. Not 
surprisingly, it is impossible to formulate a general no-hair theorem for arbitrary matter 1 . 
One needs to make certain precise assumptions about properties of matter and how this 
latter is coupled to gravity. Using these assumptions it is possible to formulate specific 
no-hair theorems only in some particular cases. For example, in the case of a general 
minimally coupled scalar field (when the Lagrangian of the scalar field is a function of 

1 When one includes a cosmological constant in the action the presence of an extra length scale makes 
no hair theorems far weaker (see for example [2]). In particular, no hair theorems generically do not work 
for asymptotically adS spacetimes. This is due to the trapping nature of adS spacetime. 
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the canonical kinetic term and the field itself), under reasonable assumptions including 
rather weak ones for the scalar potential, one can prove that scalar hair is absent for a 
static black hole, (see [3] for an overview). For a stationary black hole one has to assume 
stronger conditions for the scalar potential but even they do not cover all possibilities. 
Indeed very recently a loop-hole has been found and a hairy stationary black hole has been 
constructed numerically by Herdeiro and Radu [4, 5]. This is achieved for a (minimally 
coupled) massive complex scalar field. The crucial property is that the scalar field does 
not have the same symmetries as those imposed on space-time. This is an allowed anzatz 
for the field equations as long as the energy-momentum tensor of the scalar field does have 
the spacetime symmetries (while the scalar does not). As such a black hole is numerically 
constructed taking advantage of the superradiant instability of stationary black holes. This 
method fails for static spacetimes where a no hair theorem can be also shown to hold [6], 
[7] even though the scalar field is allowed to break spacetime symmetries. 

The situation is more complex for non-minimally coupled scalar fields, in particular, 
when the theories under consideration are allowed to be of higher order but with second 
order field equations. This is the case of Horndeski/Galileon theory [8, 9]. Due to the higher 
order nature of the theory, additional branches of solutions appear yielding sometimes to 
surprising results. Furthermore, the galileon as it is called in such theories, in general 
kinetically mixes with the graviton, therefore this opens other possibilities to evade the 
no-hair theorems. An elegant approach in order to elaborate a no hair theorem for shift- 
symmetric Galileons has been put forward in [10], where the use of shift symmetry allowed 
(again, with some assumptions) to prove such a result. 

One can, however, look for possibilities to allow for the presence of scalar hair by 
slightly twisting or bifurcating some of the hypotheses. In the case of Galileons, a non¬ 
trivial possibility in order to evade the result of Ref. [10] would be to relax the assumption 
of staticity of the Galileon field while assuming it to be true for spacetime. Indeed, for 
a standard canonical field with a potential the staticity assumption is natural, since one 
expects that at spatial infinity the scalar evolves to its minimum (as it occurs during 
cosmological expansion). In the case of a shift-symmetric Galileon, however, the static 
ansatz for the scalar field is not natural anymore, due to the non-trivial dynamics of the 
Galileon. Indeed, depending on the choice of the Lagrangian, the Galileon field may have 
several attractors for the cosmological dynamics. The non-trivial attractors — in particular, 
with de-Sitter metric — have a non-trivial evolution of the scalar field and the scalar ends 
up depending linearly on cosmic time. 

The cosmological evolution of the scalar field poses the boundary conditions for the 
galileon far from the black hole. Therefore a static ansatz for the black hole galileon 
would clearly be at odds with a cosmological galileon, since at large radii the assumption 
is clearly violated. Thus a natural choice for a shift-symmetric Galileon theory would be, 
in fact, to assume that the scalar field is time-dependent. The simplest way to implement 
this assumption is to take a linear time dependence of the scalar field 2 . This assumption 

2 Due to the non-trivial mixing of the scalar and the graviton, interesting physical phenomena emerge in 
local physics, when the scalar field is time-dependent [11]. 
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obviously bifurcates the no-go theorem of Ref. [10], giving a principal possibility to create 
a non-trivial configuration of the scalar field. 

On the other hand, one can look for a loophole in the proof of a no-go theorem. For a 
shift-symmetric Galileon, an important step of the proof relies on the physical requirement, 
that the current associated with the shift-symmetry, is non-diverging at the horizon. This 
means that the radial component of the current, J r , must be zero. The non-trivial step in 
Horndeski theories is to then prove that J r = 0 leads indeed to a trivial configuration of 
the scalar field [10]. This is not trivial and in some cases as we will see simply not true 
since the theory at hand is of higher order and other possibilities do exist. 

In Ref. [12] both loopholes have been investigated to find black hole solutions with non¬ 
trivial scalar field (see also [13] for a complimentary discussion). A simple model has been 
considered: the standard kinetic term plus a Galileon term ~ A combination 

of these terms allows for the r-component of the current to have zero value, while at the 
same time the scalar field is not trivial. Moreover, the scalar field itself has been found to 
be regular at the horizon (which is a much stronger condition than the finiteness of the 
norm of the current), for the time-dependent configurations 3 . The same strategy later has 
been applied to find a number of solutions in other Galileon-type models, see e.g. [18-20]. 

In this paper we extend the work [12] to include a gauge field. In Sec. 2 we write down 
the full action and its equations of motion. Then we assume a specific spherically symmetric 
ansatz. This ansatz is static for the metric and the gauge field, but time-dependent for 
the scalar field. In Sec. 3 the system of differential equations is completely integrated, 
leading to two algebraic equations. This system of algebraic equations fully describes the 
solution, though in an implicit form. To find explicit and interesting solutions we consider 
judicious choices of parameters and integration constants, for which we are able to write 
down the solutions analytically. In the last section we take a step back and discuss the 
generality of our scalar field Anzatz and most importantly its mathematical consistency. 
In other words: is our scalar (time-dependent) anzatz just a way to obtain particular 
very special solutions for special theories or is there a deeper and more general reasoning 
behind the integrability properties? This, we believe is a crucial and central question in the 
further study of Galileon black holes. To answer this question, we consider a general scalar 
tensor theory with translational invariance. We assume a static and spherically symmetric 
spacetime along with a linear time dependent (plus radial) scalar field anzatz. We then 
prove that the equations of motion are not overdetermined and the anzatz is in fact a 
consistent one leading always to an ’’integrable system”. Most importantly we spotlight 
an intrinsic property of the system of equations that pinpoints the J r = 0 property of 
the scalar field. The time-radial component of the metric equations of motion (the flux 
of the scalar) is actually proportional to the radial component of the Noether current J ri 
associated with the shift-symmetry. Thus setting J r = 0 (with a non trivial scalar) allows 
the immediate resolution of 2 held equations rendering the system not only consistent but 

3 A number of works studied similar theories, finding solutions with a static scalar field, [14-16]. However, 
unlike the solution of [12], the scalar field explodes at the horizon for these solutions. In [17] the authors 
constructed wormhole solutions with the same type of the Lagrangian. 
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more tractable and in certain cases completely integrable analytically. We finally discuss 
the obtained results in Sec. 6 . 


2 Field equations and Ansatz 


In this Section we shall discuss the choice of the action, the field equations, as well as the 
ansatz we will opt for. The action we consider includes the metric, one real scalar held 
and one Abelian gauge held. The latter part of the action which is new in this paper is 
constructed in two steps, described in detail below. We consider a static and spherically 
symmetric ansatz for the spacetime metric as well as the gauge held. We will crucially 
consider an additional (to the radial) linear time dependence for the scalar held which will 
be primordial for the regularity of the black hole solutions we will find. The consistency 
of this ansatz will be discussed in the hnal part of the paper, since it is true for a general 
Galileon/gauge theory admitting translational invariance. 


2.1 Construction of the action 


We consider a four-dimensional action, which depends on a metric, one scalar held and 
one 17(1) gauge held. The action we will take into account can be thought of as a two-step 
construction, following the scheme, hrst presented in [12]. Let us start from the standard 
action, containing the Einstein-Hilbert term, a cosmological constant and the standard 
Maxwell gauge held, 


So 


y/^gd 4 x 


R-2A 





( 2 . 1 ) 


where EW = d^A u — d v A^ is the held strength. The variation of each of the term in 
the above action with respect to the metric gives correspondingly (up to constants), the 
Einstein tensor the metric g^ u , and the energy-momentum tensor of the Maxwell held, 

T^ ] = \ (f ii(J F v ° - 1 g, u F a gF a ^ . ( 2 . 2 ) 


We build the scalar interaction terms in the action by contracting the result of the 
variation of ( 2 . 1 ) with the tensor made out of the derivatives of the scalar held. 

With this procedure, we obtain, 


R^G lxu V»cj)V'(t>, A —»• {d(f>) 2 , (2.3) 


Now, allowing arbitrary coefficients in front of the obtained terms, the full action which 
we study in this paper reads 


S[<7/i vi </*) Afy — 

R-2A-jF^F^ + p GVV'VVV - g {defy 2 - 7 V^V> , 

(2 - 4) 

where r/, /3, 7 and A are constants and T^ f> is given by (2.2). We omitted the overall 
Alp (where Mp is the Planck mass). Since we are only interested in classical solutions, 
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this overall constant does not enter the final result. In our conventions the scalar field 
is dimensionless, as well as the gauge field. With these definitions, 77 is a dimensionless 
constant, (3 and 7 are constants with the dimension [length] 2 and the dimension of A is 
[mass] 2 . The case (3 = 0, 7 = 0 and g = 1/2 corresponds to the canonical minimally coupled 
scalar field, and minimally coupled Maxwell electromagnetic field. The term (3 
(which we will refer to as “John” term, following [21]), gives a non-minimal kinetic coupling 
between the scalar field and gravity. The theory without the Maxwell kinetic term and 
with j3 7^ 0 and 7 = 0 has been considered in [12], where various black hole solutions were 
presented, see also Refs. [18-20] for other solutions in arbitrary dimensions with different 
asymptotics. The last term in (2.4) mixes the scalar field and the gauge field. 

Note that the action (2.4), although higher order in derivatives, gives second-order 
equations of motion, as we will explicitly see later. The presence of the scalar field in 
(2.4) modifies the dynamics of gravity. This scalar field is an ingredient for modification 
of General Relativity, like in a generic Horndeski model. The gauge field can be thought 
of as an extra modification, in which case it is unrelated to the electromagnetic field of the 
standard model. 

On the other hand, one can think of F^ u as of the standard electromagnetic field. In 
this case, one should set 7 = 0 (or 7 must be very small), otherwise photons would decay 
into the scalar field, which would be detectable by various experiments. In this context, 
the action (2.4) can be seen as a Horndeski scalar-tensor theory plus matter, where the 
electromagnetic field plays the role of matter. Note that in this case, certain electrically 
charged solutions have been reported in [16]. 

In the first part of the paper we will keep our discussion as general as possible, not spec¬ 
ifying the nature of the gauge field. In the second part we give particular solutions, which 
usually corresponds to particular choices of 7. Depending on this choice, a corresponding 
interpretation of the gauge field should be given. 


2.2 Field equations and the ansatz 

Let us explicitly write down the equations of motion for the action (2.4). The scalar field 
equation, obtained by varying the action with respect to <fi, reads, 


V„ 


G"- r] <r-'lT § [) ) V„ 


= 0 . 


(2.5) 


It is interesting to see that (2.5) can be rewritten in terms of a conserved current equation 
as 

r = 0 , ( 2 . 6 ) 

where the current is defined and given by 



This property is a consequence of the symmetry of the action (2.4) with respect to the shift 
symmetry (j) —>■ <^+const. The quantity is naturally the Noether current associated with 
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this symmetry. Variation of (2.4) with respect to A^ gives the gauge field equation, 




F> JU - 7 ( ^ F» v V CT 0 + (F a - F a v V^) ) V 1 


= 0 . 


( 2 . 8 ) 


And, finally, varying (2.4) with respect to the metric gives the modified ’’Einstein” equa¬ 
tions casted in the following form 


n — vd 1 ) i vd 2 ) 

'- jr A u/ ± iiv ' A fiu i 


(2.9) 


where we defined, 


= -a v + tW 

T$ = - 2Va 0V (m ^ } - - (V„VV)(V„V A 0) + (V^V^p^ 


+-G^(V0) 2 — 


H-7 

2 ' 


i(V A V^)(V A V^) + i(D^) 2 -V A 0V p ^ }, (2.10) 

F lu ,F Ifp V r <l>V p <j> + -±g tu/ F Pff F T a Vf i <l>V r <l> 


+ 2 rjd^dycl) - rig^dvcpd^cl). 


Note that T^J corresponds to the A-term and the standard Maxwell held, as defined in 
(2.1), while T)v comes from the the last three terms in (2.4). 

The above held equations look and are very complex. In the next subsection, we, 
however, choose an ansatz, for which the equations of motion are greatly simplihed. Many 
components of the Einstein equations are trivially zero due to the symmetries of the ansatz. 
The same is true for the vector equation (2.8). Moreover, as will be shown in a later section, 
we will establish a relation between one of the components of the Noether current with one 
of the components of the metric equations. 

We assume the following static spherically symmetric ansatz for the metric, 


dr 2 


ds 2 = — h(r ) dt 2 + —— + r 2 (d0 2 + sin 2 (0)dtp 2 ), 

f{r) 

while the scalar held is taken to be linearly time-dependent as 


( 2 . 11 ) 


< P(t,r) = qt + i’(r). 

Finally, we opt for a standard “dyonic” ansatz for the gauge held A 

A^dx^ 1 = A(r)dt + B(9)dp. 


( 2 . 12 ) 


(2.13) 


Note that as long as B{6) is a general function of 9 , the ansatz (2.13) is not spherically 
symmetric. However, there is a particular choice of B{6 ), for which A^dx^ is spherically 
symmetric. Indeed, the held strength two-form reads, F = ^ A'(r)dt A dr + ^ BQ(0)d0 A dtp, 
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where Bg(9) = dB(9)/d9. The first term of this expression only depends on the radial 
coordinate and corresponds to the electric field. The second term is spherically symmetric, 
only if Bg(9) oc sin 9. because in this case Bg(9)d9 A d(p is proportional to the element 
of the solid angle (i.e. the element of the area on a unit sphere). In the Sec. 3, we will 
start from the ansatz (2.13), and we will see explicitly that the equations of motion require 
B{9) oc cos 9 reducing the ansatz (2.13) to 

A^dx^ = A(r)dt — P cos(9)d(p, (2-14) 


where P is a constant. 

3 Solutions in an implicit form and their regularity 

In this section we solve the set of field equations (2.5), (2.8) and (2.9), assuming the 
ansatz (2.11), (2.12) and (2.13). In the first part of the Section, we give the general 
solution written in an implicit from. An auxiliary function of the radial coordinate is 
introduced to simplify the problem, and this auxiliary function turns to be a solution of 
a fifth order algebraic equation. The metric functions, the scalar field and the gauge field 
can be deduced directly from this auxiliary function. For special choices of the action 
parameters and the integration constants, it is possible to reduce the algebraic equation 
on the auxiliary function to lower order and write down the solution explicitly. In the 
phenomenologically interesting case, 7 = 0, we present a perturbative solution valid for 
small electric and magnetic charges. 

3.1 General solution in an implicit form 

We take the ansatz as in the previous section, i.e. in the form (2.11), (2.12) and (2.13). 
We will see in a moment, that the equations of motion reduce (2.13) to (2.14). Indeed, 
defining 

F(r) = A', C{9) = B e , 

the Maxwell equations (2.8) yield the following PDE, 



= —(—) • 
sin(0) Vsin (9) J g 


As we will see in Sec. 5, the (fr)-component of the Einstein equation is always satisfied if 
J r = 0, i.e., 

P G rr -V 9 rr T{m) = 0. (3.2) 

Demanding then that (3.2) is satisfied, the separability constant C in (3.1) is found to 
be zero. As a result, the vector potential takes the form (2.14), as we expected. We 
note that, by satisfying (3.2), we solve simultaneously the (tr)-component of the Einstein 
equation and the scalar field equation as we will see explicitly and in full generality in the 
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forthcoming section. In other words the Anzatz (3.2) that we impose is actually not only 
judicious but also necessary for a full solution to the field equations. At this point we can 
define our auxiliary function, setting 

P(rh{r)y + lr 2 F 2 


S(r) = 


r)r 2 + 0 — 


7 P 2 


4 r 2 


so that Eq. (3.2) reads, in terms of S(r), 


fi r ) = 


h i r ) 

S(rY 


(3.3) 


(3.4) 


Our aim now is to give all unknown functions in terms of the function S. Given that (3.2) 
guarantees simultaneously that the scalar field equation as well as the (tr) equation are 
satisfied; hence it only remains to solve the (r, r) and (t, t ) components of the metric field 
equations. The (r, r )—component is an algebraic equation on if'/ yielding 


w) = 


q 2 ^4 /3 h! (r) + 7 vF (r) 2 ^ r 3 


(4 r\ ?' 4 + 4 f3 r 2 — 7 P 2 ) h(r) 


h{r) 


(7 — P) Of (r ) 2 [ P 2 (7 - (5) - 4 r 4 (77 + /3 A)] S (r) 


(3.5) 


P 0 J 

while the (t, t)-component of the metric equation can be integrated once and written as 


0 


q 2 P-jph-P)F 2 — S(r) 


(7? _/3A) r 2 + 2 / 3__Lp2 (/3 + 7) 


+C 0 S(r ) 3 / 2 


■qr 2 + 0 — 


7 


P 2 


4r 2 


(3.6) 


= 0 , 


where Co is an integration constant. 

The only remaining equation is the time component of the vector equation (3.1), which 
reads, for our ansatz, 




i + 2l sw?-\ 


!\2 9 


= Q, 


(3.7) 


where Q is a constant, related to the electric charge. 

Using Eq. (3.6) and its first derivative, the electric equation (3.7) can be rewritten as, 


0~1 + F{r) _ Q 


0 


2 ’ 


(3.8) 


Sir) 1 / 2 

We can now treat (3.8) and (3.6) as two coupled algebraic equations with respect to S and 
F yielding a full solution to the system. 

Equivalently we can simply substitute the above expression (3.8) into (3.6) in order to 
get an algebraic equation of the fifth order in S l / 2 {r) given by 


(Co7V / MO + 2/3-2 7 )' \0 2 Q 2 


(77 — /3 A) r 2 + 2 /3 — 


F 2 (0 + i) 

4 r 2 


S(r) 


+C 0 Sirf' 2 (7 r 2 + 0 - | - (7 - 0) Q 2 Sir)(3 2 = 0. (3.9) 
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The system is therefore integrable in terms of a single algebraic equation on S and the full 
ensemble of solutions can be obtained, which will depend on three integration constants 
Co, Q and P. The integration constant Co determines the asymptotic behavior of the 
solution. Indeed note that (3.6) is independent of the mass of the black hole solution since 
S(r) depends on the derivative of h(r). The integration constant Q is related to the electric 
charge, while P is the magnetic charge of a black hole. We have therefore reduced the full 
field equations to the solution of a single algebraic equation of the 5th order with respect 
to the auxiliary function S. 


3.2 Regularity for the implicit solutions 

We now argue, that a solution, given by the equations (3.6) and (3.8), generically contains 
a branch, which is regular at the future black hole horizon. We will now show that if 
r = t'h is an event horizon with h{rn) = 0 then the solution will be regular if S(r) > 0 for 
all r > rn- 

Indeed, following [12], we introduce the ingoing Eddington-Finkelstein coordinates, 

v = t + J (fh)~ 1/2 dr. (3.10) 

In the coordinates (v,r), the metric takes the form, 

ds 2 = — hdv 2 + 2 y/h/f dv dr + r 2 dfl 2 . (3-11) 


The ratio h/f is simply S from (3.4), implying that the metric is regular if S(r ) fulfills 
the above property. Similarly, from (3.8) one can see that F(r) is regular, unless the 
expression in the parentheses is zero at the horizon, which, however, requires some tuning 
of the parameters of the model. 

One may also worry about the regularity of the norm of the scalar current, J 2 = 
Since the action is shift symmetric, the norm of the current is a natural quantity 
to check. For our solution, the radial component of the current is zero by construction, 
see Eq. (3.2). Thus, the only non-zero component of J^ comes from the time component. 
From the definition of the current, Eq. (2.7) and using (2.12), we have, 


J 1 = 



(3.12) 


Note the presence of g u = —l/h(r) in the above expression, which is clearly divergent 
at the horizon and may be dangerous for regularity. However, other terms also contain 
singular parts, which combine to a non-singular expression. Indeed, calculating explicitly 
G u for the ansatz (2.11) and using (2.2), one finds, 


J t 


1 

r 2 h 2 


h 




PKfr)' 


XT r2 fA ' 2 


q, 


(3.13) 


which, by virtue of (3.3), (3.4) and (3.13) becomes, 

rt Pq S'(r) 
r S(r) 2 


(3.14) 
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Again if S(r) is smooth and non-zero at a horizon then the norm on the current is not only 
finite, but it is in fact equal to zero, 


lim J 2 = 0. (3.15) 

r^rr H 

It is interesting that not only the norm of the Noether current is regular, but also the 
scalar field itself again under the same condition on S. Let us see how this works out. Near 
the horizon, r —>• m, h{rjj) —>• 0, and using (3.3), one can rewrite (3.5) as follows, 

<=r H = ±{s{r) 1 / 2 + 0(\). (3.16) 

This is of course singular and one would immediately conclude, if the scalar did not have 
a time dependance, that there is a scalar singularity there. However, substituting this 
expression (with the plus sign) into (2.12) and taking into account (3.10), one finds, 

J dr +0(1). (3.17) 

The expression inside the integral of (3.17), naively, would give a divergence of <j) at the 
horizon, since it contains / and h in the denominator. However, using (3.4) one can verify 
that the two terms inside the integral cancel each other, thus yielding 4> = qv + 0(1) at 
the horizon. Thus the scalar held is not diverging at the future black hole horizon! 

However, the scalar held, although regular at the future horizon of the black hole, 
usually diverges at the de-Sitter horizon, provided that the solution is asymptotically de- 
Sitter. This happens because the coordinate v is not regular at the future de-Sitter horizon 
(independently of mass). Instead, the outgoing coordinate u = t — f (//i)~ 1//2 dr is regular 
there 4 . The solution (3.16) with minus gives a regular solution for cj) at the future de- 
Sitter horizon. We should however warn the reader that the action itself depends only on 
derivatives of the scalar field and not the scalar field itself and for the derivative field all 
quantities are finite and well-defined. 

4 Particular explicit solutions 

The master algebraic equation (3.6) is of fifth order, and writing down a full solution 
involves complicated integrations which means that one can not write down the general 
solution explicitly. Nevertheless, for some special cases the solutions can be written down 
explicitly, as we will now see. 

4.1 Case Co = A = rj = 0 

It is easy to note that the equations (3.6) and (3.8) are considerably simplified for Co = 0, 
since in this case (3.8) implies that F(r) oc 5(r) 1//2 , and substituting this into (3.6), the 
latter becomes a linear algebraic equation on S(r). As a direct consequence S(r) is found, 
and it is not difficult to write down the solutions for all other functions. However, the 

4 We thank T. Kobayashi and N. Tanahashi for pointing this out to us. 
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full expressions are still somewhat cumbersome; therefore, we present here the solution for 

A = r] = 0. 

The metric functions as well as scalar field profile take the following form 


,,, , H "v/2 [P 2 (7 — /?) — (/? — 3 7 ) Q 2 ] ^ 

h(r) = 1 - ; _ - arctan 

r 4 r y//3 [Q 2 (7 — /3) — P 2 (/3 + 7 )] 


2y/2pr 


VQ 2 {i-P)-p 2 (18 + 7 ) 


f(r) = 1 + (7 - /3) 


Q 2 P 2 h + P) 


h(r), 


(^'( r )) 2 = p r 2 2 f{r) 2 (W ~ /( r )) + 


8 r 2 P 8 /3r 2 

Q 2 {l ~ P) p 2 (7 + P) 


Pr 2 f{r ) 2 
together with 

F(r) = ±-|(l + 


8 /? 


8/3 


(4.1) 


0 2 (7-/3) P 2 (/3 + 7 )V 1/2 

8r 2 ^ 8/3 r 2 

aled electric cha 
we also fixed q to have a particular value, 


, Fev = C{6) = P sin(0), (4.2) 


where we defined the rescaled electric charge Q as, Q = A In the above expressions 


g 2 /3 = 2, 


(4.3) 


such that h{r) —>• 1 for asymptotically flat solutions, so the time coordinate is defined 
properly and we don’t have a solid deficit angle. In the above expressions we also implicitly 
assume that 

Q 2 (l ~ P) > P 2 {P + l). (4.4) 

In this case, the solution (4.1) and (4.2) corresponds to an asymptotically flat black hole 
with S strictly positive so that the scalar field is well defined beyond the event horizon of 
the black hole given by the condition h{rn) = 0. Both the metric and the electromagnetic 
field are modified as compared to the charged black hole solution in GR. Far from the black 
hole, that is for r —>• 00 , the electric held (4.2) has its usual form, F ~ 1/r 2 , however, with 
an effective charge Q. On the other hand, the magnetic held is not modihed. It is also 
worthwhile to note that the mass of the black hole as measured at infinity is not given 
simply by fi, but it is modihed, since the last term in the hrst equation of (4.1) behaves as 
1/r at r —>• 00 . 

One can also check that the norm of the Noether current J 2 and the gauge held, (4.2) 
are regular at the horizon. The scalar held (4.1) (for the branch with the plus sign) is 
regular at the future black hole horizon. 

When the condition (4.4) is not met, i.e. for Q 2 (7 — /3) < P 2 (/3 + 7 ), Eq. (4.2) 
becomes imaginary at small r, thus the solution is not physically viable and we do not 
discuss it. 

By setting P = 0, we easily recover the purely electric case, with no magnetic charge. 
The condition for a well-defined asymptotically hat solution is then simply 7 > /3. The 
purely magnetic case is also easily obtained from the above solution by setting Q = 0. The 
parameter 7 in this case must have a negative value, 7 < —/3. 
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Let us now consider another simple case for which the two metric functions are equal 
as it occurs for the standard Reissner-Nordstrom solution in GR. 

4.2 Case f(r) = h(r ) 

The form of the expression for S(r), see Eq. (3.3) suggests another explicitly solvable case. 
Indeed, assuming the two metric functions to be identical, f(r) = h(r ), then we have 

S(r) = 1. (4.5) 


Using (4.5), one easily finds that 

h(r) = 1 — — + 


fi , ?jr 2 7 ( Q 2 + P 2 ) 


3/3 


+ 


4/3 r 2 


, ,// ^2 1 - f( r ) 2 

( r )) = t ,, 2 q ’ 


fir) 


Q 


F tr = F(r) = Fq = C {9) = P sin (9). 


(4.6) 


The coupling constants, the constants of integration and the the “velocity” q are related 
as 


P 2 P (A 7 + ry) = Q 2 rj (7 — /3), q 2 = C 0 = - ( V - A). (4.7) 

PV V 


The electromagnetic field is the same as in the case of the Reissner-Nordstrom solution. 
The metric of this solution is identical to the dyonic version of the Reissner-Nordstrom- 
de-Sitter metric, however, up to the redefinitions of Q and P. One can verify that the 
metric, the gauge field and the norm of the Noether current are regular at the horizon. 
The scalar field is regular at the future black hole horizon for the branch with the plus 
sign, but diverges at the future de-Sitter horizon. 

When the magnetic charge is zero, P = 0, the first equation of (4.7) implies /3 = 7 . 
On the other hand, when the electric charge is zero, Q = 0, then the same equation leads 
to the relation A 7 + 7 = 0, but in this case, the time-dependent character of the scalar 
field is lost. 

This solution generalizes in a rather elegant way the self tuning black hole [12] in 
the presence of an electric or magnetic charge. The self-tuning of the bulk cosmological 
constant is robust not only at the presence of gravitational mass but also in the presence 
of electro/nragnetic charge. 

Having found this stealth Reissner-Nordstrom-de-Sitter metric we can now find with 
relevant ease new solutions by a simple observation [22]. Let us assume for simplicity that 
P = 0 in which case this class of solutions dictates that (3 = 7 . Then we have immediately 
that the electric strength has the usual Maxwell profile, F(r) = % and the equation for 
S (3.6) is a third order polynomial equation in \/S. Given that y/S — 1 is a root of this 
polynomial we simply divide the two, obtaining a second order polynomial equation for 


Vs, 


s- 


u 


R 2 + 1 


s 1 / 2 — 


u 


R 2 + 1 


= 0 


(4.8) 
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where 


r] + f3 A 
rj — (3 A 

and we have rescaled the radial coordinate to R = \J^ r - It is then a straightforward 
exercise to read off the charged black hole solutions corresponding to the two remaining 
roots, 


h(R) = u— 


M yQ 
R + 4f3R 2 ' 2 R 


2 u 2 

+ — 


arctan R ± arctan 


Ry/u 


V4 R 2 + 4 + 


with f(R) = where S is the square of the solution(s) of (4.8). The scalar field can 
be directly calculated from (3.5) and is regular. This solution is the charged version of the 
black hole solution found in [ 22 ], Note that for A = 0 we have u = 1, and we have well 
defined asymptotics for h. Note that for large R we have, 


h(R) ~ 1 - 


C n In R, 
R 




with Ci = Ci(u, n, Q)\ in other words we have no longer a standard Newtonian fall-off and 
the mass and charge terms are rescaled. Also we note that f(R) ~ R 2 , and hence our 
spacetime asymptotes an Einstein static universe. 


4.3 Perturbative solution for 7 = 0 and 7 = A = 0. 

In this subsection we consider the phenomenologically interesting case 7 = 0 where the 
electric field is a Maxwell field which is not coupled to the scalar Galileon field. We 
could not find simple enough explicit solutions to Eqs. (3.6), (3.8). Instead, we perform a 
perturbative analysis for small charges Q and P. For P = Q = 0 and 7 = 0 , Eq. (3.8) is 
satisfied for F = 0, that is the gauge field is identically zero. Eq. (3.6) is solved then for 
S = So, where So is a constant, which depends on Co, /3 and q. Substituting S = So into 
(3.3) we find, ho(r) = So — where g is a constant of integration. To have canonically 
defined time at spatial infinity we have to set So = 1, which fixes Co from (3.6) as follows, 
Co = 2 — q 2 /3. From (3.4) one can easily see that 

fo{r) = h 0 (r) = 1 -^-, (4.10) 

and from (3.5) we have, 

$) = ± J~\f'K- ( 4 - n ) 

ho v 

This solution is nothing but the stealth solution defined on the Schwarzschild metric found 
in [12]. Now we assume that the charges Q and P are non-zero, but ’’small”, and we find a 
perturbative expansion in Q and P around the solution (4.10), (4.11). Eq. (3.8), after the 
redefinition Q —>• Q(3 2 , implies the standard Maxwell law, 

F = % < 4 - 12 ) 
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Substituting (4.12) into (3.6) one gets a correction to S(r), 


•SM = 1 - ( 2^2 


1 \ Q 2 + P 2 


3 q 2 /3 - 2 ' 


(4.13) 


linear in Q 2 + P 2 the strength of the EM field. It is not difficult then to find the leading 
order corrections to the metric, 


u Q 2 + P 2 l 

r + 2(3q 2 /3 — 2) r 2 ’ 


/ = M 7 ") 



Q 2 + P 2 \ 
2r 2 (3q 2 {3 — 2) / ’ 


(4.14) 


and to the scalar field, 


if) 1 = ± 


q\frh! 



Q 2 + P 2 \ 
8f3q 2 r 3 h') ' 


(4.15) 


The metric (4.14) is asymptotically flat. It is not of the Reissner-Nordstrom form, since 
h(r ) ^ /(r) at the leading order in Q and P. The solution (4.12), (4.14) and (4.15) is valid 
as long as the charges Q and P satisfy the relation, 


Q 2 + P 2 <2r 2 |(3g 2 /3 —2)|. 


(4.16) 


which is true generically before we hit the event horizon of the black hole. Again, one can 
easily check that the solution (4.12), (4.14) and (4.15) with the plus sign is regular at the 
future black hole horizon. 


4.4 Static case for f5 = 7 


Let us now consider the case q = 0. The scalar field then has no time dependence and as a 
result we expect for the scalar field to be singular at the horizon. However, these solutions 
are not immediately rejected because the action depends only scalar derivatives and they 
are in fact regular. So we study briefly this case here. Previously, progress in this direction 
was reported in [16] where the authors considered the static q = 0 case but for a minimally 
coupled gauge field 7 = 0 . Here we will consider the case /3 = 7 . The auxiliary function 
S(r) takes the form, 


4 [2 (A p - 7 ) r 4 — 4 p r 2 + P 2 /3 \ 2 
C 0 2 (4 77 r 4 + 4 /3 r 2 — P 2 (3) 2 


For simplicity we will consider only the electric case and set P = 0 before evaluating the 
relevant functions. The dyonic or magnetic solution are obtained with ease using the above 
5. 


For P = 0 we have, 


h(r) 

f( r ) 

W(r )) 2 


C ri-pA.fr 2 t (3 7 + /3 A) (77 — /3 A) i /? (77 + /3 A ) 2 _^_ 

3 7 Co Co 1 1 2 

Cj±V^+0f w , 

[(77 — /3 A) r 2 + 2 f3} 2 /3 
r 2 (77 + /3 A) 

/3 (rj r 2 + j3) f(r)’ 




? 

r 


(4.18) 
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Note here that the scalar function t/>(r) will be singular since there is no time dependence 
to cancel out the radial singularity. However the norm (V0 ) 2 appearing in the action is 
regular and finite at the horizon. The electric held strength is given by, 

F rt = B[r) = (4.19) 

and it is noteworthy that the held is identical to a Maxwell held for spherical asymmetry 
unlike what happens in the 7 = 0 case. 


5 Relation between the (tr)-component of the Einstein equation and the 

Noether current. 

As we mentioned in the introduction we would now like to take a step back and prove 
that our ansatz (2.11), (2.12) and (2.14) is mathematically consistent. This is not trivial 
since the scalar held does not obey the same symmetries that the spacetime. It turns out 
that there is an interesting relation between the (fr)-component of the metric equation 
associated with the hux of energy of the scalar held and the Noether current associated 
with the shift symmetry of the action. In this section our hypothesis is the above ansatz 
but we do not require a specihc form of the action (2.4), but rather consider a general 
action given by 

5 = y'd 4 x£(^,V^,V 2 0,F /il/ ), (5.1) 

In the above we assume that the action does not depend on cj) explicitly, but only on its 
derivatives (as is the case for (2.4)); this requirement ensures the existence of a Noether 
conserved current central to our whole discussion 5 . We also require that the vector held 
is of Maxwell type, i.e. it enters only via the gauge invariant combination = d^A v — 
d v A^. These properties are crucial for finding black hole solutions with non-trivial scalar 
conhgurations as we saw earlier on. Let us analyse the structure of the equations of motion 
for the ansatz ( 2 . 11 ), ( 2 . 12 ). 

The equations of motion for the action (5.1) can be written generically as, 


gW = 


g(a) 


1 5C 

1 5C 

V~9 SA v 

_2_ 5_C_ 

V~9 Sg^ 


= 0 , 


= 0 , 


= 0 . 


(5.2) 

(5.3) 

(5.4) 


For the ansatz (2.11), (2.12) and (2.14) the non-trivial components of the metric 
equation are the (ft), (tr), (rr) and [99) components (the [tp<p) component is also non¬ 
trivial, but it is a simple consequence of the [99) component, due to spherical symmetry, 
£q 9 q sin 2 9 = £<^). Thus, there are four metric equations. There is one scalar equation, 
£{<t>) = 0 and one non-trivial component of the gauge held equation, the time component. 


5 The expression for Noether current in the Horndeski theory can be found in [23]. 
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Thus, there are in total six equations. Now notice that the ansatz (2.11), (2.12) and 
(2.14) contains only four independent functions, namely, /(r), h(r), 4>{r ) and A(r). This 
a priori means that the system of equations is overdetermined. However, as we will now 
demonstrate, not all the equations are independent. 

Let us consider the infinitesimal coordinate transformation, 


—>■ x^ = (5.5) 

The scalar field in the new coordinates is 4>(x) = (j>(x) which, expanded in terms of x^, 
gives, 4>(x) = 4>{x) + cj)^ + £>(£ 2 ), yielding 

4>(x) = <j>{x) - (5.6) 

to leading order. Similarly, for the metric we have, g ,JjV {x) = g IMl '(x ) + g ail £, v a + g au £%, and, 
on the other hand, g^ u {x) = g^(x) + dgg IM '^. which combine to, 

r u (x) = !T(x) + = sT{x) + £* v + ^. (5.7) 

Under the infinitesimal change of coordinates (5.5) the vector field transforms as, 

\(x) = A m (x) - (£ a A a ) iAt . (5.8) 


Since the action (5.1) is invariant under the coordinate changes, we obtain, 

SS = ~J JxV^gSM KC- Jd A x^ i £^ A) (£, a A 0 ). IJ + Jd*xV=gS$e i, ' = 0 . 

(5.9) 


Integrating by parts one can get rid of the second term since 

V m £ ( % = -J=dt{V=gSl A) ) = 0. (5.10) 

Taking this into account and integrating by parts the last integral (5.9), we obtain, 


= o. 


(5.11) 


for arbitrary Now, Eq. (5.11) consists of two non-trivial relations, the time and radial 
components, 

rj£< * ) + ~7^3 dr ( v/= ® 9 " £ '<” ) ) = °’ ( 5 - 12 ) 

W + ^=a„ (v 9 =»9'“'4? ) ) - = °- < 5 - 13 ) 

As we can see, two among the six non-trivial equations of motion are not independent. 
E.g., using (5.12), one can express the scalar equation in terms of the (fr)-component 
of the Einstein equation. Then, the (00)-component can be found from (5.13) in terms of 
(fr), (tt) and (rr) components of the Einstein equation. We have four independent ordinary 
differential equations with four functions, hence the system of equations is mathematically 
consistent. 
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This is sufficient but we can obtain an integral version of (5.12), which nicely connects 
the Noether current associated with the shift symmetry (j> —>• </>+const with the scalar 
flux £$. Indeed, since the action (5.1) does not explicitly depend on (j), but only on its 
derivatives, the equation of motion for the scalar field can be written as = 0, 

where is the Noether current. Therefore, substituting 8^ = into (5.12) 

and integrating over the radial coordinate, one finds, 

qj r + g rr S[f = const. (5-14) 

The constant in (5.14) is yet to be determined. It is not arbitrary and in fact, a more 
careful investigation shows that it is zero. Indeed, from (5.9) we obtain, integrating by 
parts (and taking into account (5.10)), 

I d A x^g + I ^xy/=gSjS)^ = 0. (5.15) 

To flesh out the result let us consider a particular infinitesimal change of coordinates, 
see (5.5), 

^ = (£(r), 0,0,0). 

Taking into account that for the assumptions above the only non-zero component of v 
are £ t,r and £ r,t 1 Eq. (5.15) can be rewritten as, 

J d A xyf=g J r d r (g£(r)) + j tfxy/^gE^ (£* ;r (r) + = 0. (5.16) 

Now, remarkably, for the ansatz above, one finds £ t;r (r) + £ r; *(r) = g rr £, r (r). Thus from 
the last equation we obtain, 

- q.r = 8l 3) g rr . (5.17) 

The above equation (5.17) is identical to (5.14) for const=0. We have thus demonstrated 
that our time dependent scalar field Anzatz is in fact compatible with the absence of flux 
for the scalar field and thus a general static and spherically symmetric metric. 

6 Discussion 

In this paper we discovered new exact black hole solutions in a scalar-tensor-vector theory, 
given by (2.4). The action, besides the Einstein-Hilbert term, consists of the standard 
kinetic term for the scalar field, the standard kinetic Maxwell term for the gauge field, 
and also two non-linear terms, describing a kinetic (non-minimal) mixing of the scalar field 
with the graviton, and a kinetic mixing of the scalar field with the gauge field. The action 
is shift-symmetric in the scalar held, i.e. the action depends on the derivatives of the scalar 
held, and not on <f> itself. Also, the vector held enters the action only through the Maxwell 
combination — dyA^), therefore the action is gauge invariant. 

We searched solutions for a static and spherically symmetric ansatz for the metric 
(2.11) and the gauge held (2.14), while for the scalar held we crucially assumed a time 
(and radial) dependent ansatz of the form (2.12). We showed that this ansatz is consistent 
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as long as the time dependence of the scalar is linear. The equations of motion (2.5), (2.8) 
and (2.9) reduce to four independent ODEs with four functions entering the ansatz. We 
also we were able to prove that for this particular ansatz the r-component of the Noether 
current, associated with the shift symmetry, is proportional to the (fr)-component of the 
Einstein equations, see Eq. (5.17). 

The system of the field equations for the precise ansatz we consider turns out to be 
integrable: we obtained two algebraic equations (3.6) and (3.8), which fully determine the 
solution. Having found the auxiliary function S(r) and the electric held, F(r) from this 
algebraic equations, other functions of the ansatz can be calculated from (3.3), (3.4) and 
(3.5). The solution contains three integration constants: Q is proportional to the electric 
charge of the black hole, P is related to the magnetic charge of the monopole, and the 
constant Co which determines the asymptotic behavior of the solution. Note that there is 
a fourth integration constant (related to the mass of the black hole), which is not present 
in the master equations (3.6) and (3.8), however, it appear when one integrates (3.3) for 
a given S(r ) to find h(r). Depending on the choice of parameters and the integration 
constants, the solution describes a black hole with various asymptotic behaviors at infinity. 
An important point is that for a generic choice of parameters the metric, the gauge held 
and the norm of the scalar current are regular at the horizon(s). As for the scalar held 
itself, its solution normally contains a branch, which is regular at the future horizon. If the 
metric is asymptotically de-Sitter, the scalar held is regular either i) at the future black 
hole horizon and at the past de-Sitter horizon, or ii) at the future de-Sitter horizon and 
the past black hole horizon. 

For some choices of the parameters of the action and the constants of integration we 
gave explicit solutions, see Sec. 3. In particular, we studied a solution for 7 = 0 . This case 
is interesting from the phenomenological point of view, since it corresponds to a minimal 
coupling of the gauge held to the gravity, therefore it makes sense to think of the gauge 
held as the electromagnetic held of the Standard model. The action (2.4) in such a case 
describes Horndeski theory with a matter held. 

Here, we have presented a quite detailed scan of the analytic solutions and it will 
be interesting to explore the thermodynamics of the configurations found in this paper 
and to compare them with the standard Reissner-Nordstrom solution for example (which 
corresponds to the case /3 = rj = 7 = 0). This issue is rendered non trivial by the fact that 
the scalar held is time-dependent, and hence as a hrst step in this task, one can analyze 
the purely static case, see Section (4.4). 

One can also explore if there exists more general Horndeski sectors with a shift sym¬ 
metry that may be integrated as in the present case. For spherical symmetry for example 
it is clear that our result in the last section is a very important step in this direction. To 
simplify the discussion let us assume no gauge held but a general shift symmetric theory. 
For the anzatz of the scalar held that we consider taking J r = 0 renders the system in fact 
integrable in the following sense. Indeed one is left with the (r, r) and (t, t ) equations to 
solve with respect to the radial scalar ip'(r) and the metric function h{r). The (r, r) is an 
algebraic equation with respect to iftir) and hence can be in principle solved algebraically. 
Then, at the end, we are left with a non-linear ODE with respect to h(r) which can be 
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solved analytically or numerically. 

Another important aspect is to study cases of lesser symmetry, for example not static 
but rather stationary solutions. There, one would like to go further and find an argument, 
such as the one for the radial part of the current being zero, which again will render the 
field equations integrable. We hope to report progress on these topics in the near future. 
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